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In a number of previous studies, we have investigated the use of the volume element of the Bures 
(minimal monotone) metric — identically, one-fourth of the statistical distinguishability (SD) metric 
— as a natural measure over the n 2 — 1-dimensional convex set of n x n density matrices. This has 
led us for the cases n = 4 and 6 to estimates of the prior (Bures/SD) probabilities that qubit-qubit 
and qubit-qutrit pairs are separable. Here, we extend this work from such bipartite systems to the 
tripartite "laboratory" quantum systems possessing U <S> U <8> U symmetry recently constructed by 
■ Eggeling and Werner (Phys. Rev. A 63 [2001], 042324). We derive the associated SD metric tensors 

for the three-qubit and three-qutrit cases, and then obtain estimates of the various related Bures/SD 
£SJ ■ probabilities using Monte Carlo methods. 
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. !• INTRODUCTION 

in ; 

In a number of previous studies [1-4] — developing upon an initial idea of Zyzckowski, Horodecki, Sanpera and 
Lewenstein [5] — we have investigated the use of the volume element of the Bures metric [6-8] as a natural measure 
£f) . on the quantum states. In particular, we have been interested in its use to evaluate the a priori probability that a 
(bipartite) state is classical in nature (separable) or, equivalently, the complementary probability that it is quantum 
(-h , (nonseparable/entangled). The Bures metric serves as the minimal monotone metric [9]. The probabilities of sepa- 
fS |' rability its volume element provide, appear to be larger than for any other member of the nondenumerable class of 
! 1 * . monotone metrics. (Classically, there is a unique monotone metric — the Fisher information metric [10,11].) In this 
sense, it also constitutes an upper bound on other (arguably plausible) prior probabilities of separability — while the 
maximal monotone metric [12,13] would provide a lower bound. (The Bures metric is associated with the symmetric 
logarithmic derivative and the maximal monotone metric with the right logarithmic derivative [14].) 

In [2], we found for certain low-dimensional scenarios involving two qubits, analytically exact Bures probabilities 
of separability. This further led us in [3] to use numerical methods (quasi-Monte Carlo integration), employing a 
recently-developed parameterization of the 4x4 density matrices based on SU (4) Euler angles [15], to investigate the 
conjecture that the Bures probability of separability of two arbitrarily coupled qubits — occupying a 15-dimensional 
convex set — might also be exact. In [16], we found numerical support for the conjecture that the Bures probability of 
separability of two qubits is 168 °^ Aa ~ 0.0733389, where CAg is the silver mean, \/2 — 1. (Our primary conjecture there 
was that the volume — as measured in terms of (four times) the Bures metric, that is, the statistical distinguishability 
(SD) metric [8] — of the separable two-qubit states is — j 2 -. Certain other monotone metrics of interest also appear to 
give volumes of the separable states that are simple multiples of crAg, such as lOcr^g for the Kubo-Mori metric [16]. 
These analyses, however, cast some doubt as to the naturalness of our earlier conjectures in [4] in regard to the Bures 
probability of separability of a qubit and qutrit.) 

The development by Eggeling and Werner [17] of new ("laboratory") systems of tripartite states — composed 
of three subsystems of equal but arbitrary finite Hilbert space dimension (d) — provides us with an opportunity 
to take this line of analysis to a finer level of description than previously, that is, by attaching Bures volumes to 
the sets of biseparable, triseparable states, as well as those having positive partial transpose with respect to the 
biseparable partition (cf. [18]). These states were devised by Eggeling and Werner to obtain a dimension-independent 
characterization of the separability properties of symmetric states. They possess an explicit parameterization as linear 
combinations of permutation operators, which is helpful in explicit computations. Also, there is a "twirl" operation, 
which brings an arbitrary tripartite state to this special subset. 

It has proved convenient in previous analyses to speak in terms of the SD metric rather than the Bures one, and 
we will follow such a course here. All the probabilities given below will, of course, be the same using either the Bures 
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or SD metric. As a point of reference, for the bipartite (one-dimensional laboratory) Werner states [19] composed of 
two qubits, the Bures/SD probability of separability was found to be exactly \ [2, sec. II. A. 7]. 

II. ANALYSES 

Five parameters (r_, r+, ri, T2, r^) characterize the Eggeling- Werner (EW) states, in general, but only four in 
the three-qubit case (where the parameter r_ degenerates to identically zero). (Also, it is convenient to employ a 
"dummy" bound variable, r = 1 — r_ — r+, as in [17], as well as utilize the transformation to spherical coordinates, 
n = R cos #,r2 = R sin cos <fi, r^ = Rsin9 sin</>.) For an EW-dcnsity matrix, the parameters obey the relations [17, 
eq. (6)], 

r+,r-,r > 0, r++r_+r = l, r\+rl+rl<rl. (1) 
The eigenvalues of an EW-density matrix arc of the form, 

r + ij.-_T • \ r - \ 1 

where 



— (multiplicity — (m. vJ), — (r ± \Jr\ + r\ + rg) (m. j/ ), (2) 



d 3 + 3d 2 + 2d d 3 -M 2 + 2d d 3 -d 
v + = g , = g , «* = —• (3) 

(The multiplicity v + is equal to the dimension of the corresponding Bose subspace divided by the dimension — that is, 
1 - of the respective [trivial] irrep, while v_ similarly corresponds to the Fermi subspace and v$ to the para-subspace.) 
So, for d = 2 there are, in general, three distinct eigenvalues, and for d > 2, four distinct eigenvalues. 

Since we were able rather readily to determine, in addition, the corresponding eigenvectors of the EW states (p) 
we examined, we could directly implement the basic formula for the Bures metric, 

i , , \2 1 V- I < Oi\& P \(3 > | 2 fA . 

a.fj ^ 

(where \a >,a — l,...,n, are eigenvectors of p with eigenvalues X a ) without having to rely upon any of the interesting 
(indirect) methods Dittmann developed [20] — and we have applied elsewhere [2,4,21] — in order to avoid the (often 
problematical) computation of eigenvalues and eigenvectors of p. (For the SD metric we replace the coefficient \ in 
(4) with 2.) We also observe the relationship 



rfs(pi,p 2 ) 2 = 2-2 v /F(p 1 ,p 2 ), (5) 

where F{p\,p2) is the frequently-used Bures fidelity between density matrices p\ and p2- [6,22,23].) 

We first examined the EW-states in the qubit (d = 2) and qutrit (d = 3) cases. It turned out that the Bures metric 
tensor in the qubit (four-parameter [r„ =0]) case is identical to the corresponding 4x4 submatrix of the Bures metric 
tensor in the qutrit (five-parameter) case. Also, numerical evidence we have adduced strongly indicates that the 5x5 
metric tensor (and hence the volume element) in the case of three four-level (d — 4) quantum systems is identical 
to that in the three-qutrit case. Presumably, this holds true for tripartite EW-states with still higher-dimensional 
subsystems (d > 4). (In fact, T. Eggeling has been able to independently confirm this proposition, as well as the 
formulas (6)-(15) given immediately below by deriving a general formula — in terms of the subsystem dimension d 
— for the Bures metric tensor of the EW tripartite states, on the basis of representation theory (App. I). He has also 
been able to show that the scalar curvature of the Bures metric is equal to 20 + 18/r , and thus diverges for m = 
(cf. [24]).) "Surprisingly, it turns out that the separability sets we investigate are also independent of dimension" [17]. 

In terms of the spherical coordinates, the SD metric tensor elements (gij) which are not identically zero are the 
simple functions, 

1/2 1 1 \ 

9r - r - = 2{- + WTR-^TR)' (6) 

9r - r+ ~ (-r -R)(-ro + RY (0 
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R . . 

(r + R)(-r + R) 

1/21 1 \ 

9r+r+ ~ 2\^ + WTR~ -r Q + RJ' (9) 



R_ 

(r + i?)Fr + J?) ' 



gr + fl = - 7: r^vT~: 5T ' ( 10 ) 



(-r - #)(-r + J?) 



ro ' 

R 2 sin 2 ( 
ro 



(12) 



(13) 



The SD volume element (Iffy I 1 ' 2 ) in the qubit case (first, having set r_ = in the above and deleted the row 
and column of the 5x5 metric tensor corresponding to r_ before computing |gij| 5 ) is (reverting to the original EW 
coordinates), 

Iffy I qubit I ; 



r v\j r +( - r i - r l - ( r o + r 3 )(-r + r 3 )) 



and, quite similarly, in the qutrit case, 



9ij\qitrit — ; 



rn\jr-r + (^- r\-r\- (r + r 3 )(-r + r 3 )) 



To normalize (14) to a probability distibution over the EW qubit-qubit-qubit states, one must divide it by 27r 2 /3, 
while to normalize (15) over the EW qutrit-qutrit-qutrit states, one divides it by 7r 3 /2. 

Eggeling and Werner also gave explicit parameter ranges, so that p : (1) is biseparable with respect to the partition 
1 1 23 ; (2) is triseparable; and (3) has a partial transpose that is positive with respect to the first tensor factor. We 
have used these explicit ranges in numerical (Monte Carlo) simulations to estimate the relative probabilities that an 
EW-state has any of these three special properties. In the qubit (r_ = 0) case, however, biseparability simply implies 
that the corresponding partial transpose is positive [17]. 



A. Probabilities for qutrit case 



For the three- qutrit case, based on some 131 million randomly generated points, an estimate of .0963689 was gotten 
for the SD/Bures probability of a positive partial transpose, .0694443 for the SD/Bures probability of biseparability, 
and .0165952 for triseparability. The sample was compiled in five roughly equal subsample. (Actually, for each 
subsamplc, one billion random uniformly distributed points — r_ and r + being drawn from [0,1] and ri,r 2 ,r 3 from 
[-1,1] — were tested to ascertain if they fullfilled the requirements (1) for an EW-state. Approximately 2.6 percent 
did in each subsample.) In Table I, we also show the results for each subsample, so one might gauge the stability of 
the overall results. (The probabilities were computed by taking the ratio of the SD volume assigned to the particular 
subset of states in question to the SD volume assigned to all those approximately 131 million states meeting the 
requirements for a general EW-state.) The standard deviations over the five samples are also given — assuming the 
nearly equal sample sizes are, in fact, strictly equal. (Since the true probabilities arc unknown, but only estimated, 

we use the bias-adjusted standard deviation, J s i=i(^'~ x ) ; w \th n = 4 rather than n = 5.) 
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Based on the results in Table I, it would perhaps appear to be somewhat of an overstatement "that the two sets 
[biseparable states and states with positive partial transposes] come to be remarkably close (see. Fig. 11)" [17, p. 
10]. 



sample size 


Bures/SD prob. of positive part, transp. 


Bures/SD prob. of biseparable 


Bures/SD prob. of triseparable 


26,176,260 


.0960822 


.0692970 


.0141430 


26,177,694 


.0962687 


.0691823 


.0142139 


26,179,742 


.0964395 


.0696499 


.0143551 


26,183,533 


.0963464 


.0693287 


.0141516 


26,173,004 


.0967076 


.0697630 


.0144027 


130,890,233 


.0963689 


.0694443 


.0142526 


std. dev. 


.00023 


.000249 


.000119 



TABLE I. Estimated Bures/SD probabilities, based on Monte Carlo simulations, that a three-qutrit EW-state: (1) has a 
positive partial transpose with respect to the partition 1|23; (2) is biseparable with respect to the partition; (3) is triseparable; 
and (4) is biseparable with respect to the partition given that it is triseparable. For each of the five subsamples, one billion 
points were initially generated, from which those corresponding to EW-states were selected. 



Zyczkowski, Horodecki, Sanpera and Lewenstein, in the conclusion of their pioneering paper [5], had asked the 
question: "Has the set of separable states really a volume strictly smaller than the volume of the set of states with a 
positive partial transpose?" So, their question is answered affirmatively in the context here. 

In our initial attempts to calculate the entries of Table I, however, we had been quite perplexed to find that roughly 
14 percent of the probability mass assigned to the triseparable states appeared not to meet the EW-criteria for 
biseparability. (A randomly generated example of one such state had parameters r+ = .27, r_ = .1, n — .589304, r2 = 
.08100014, r3 = —.138433.) After correspondence with T. Eggeling regarding this clearly paradoxical situation, he 
concluded that it was necessary to incorporate — which we proceeded to do, as reflected in Table I exhibited here — 
the additional constraint, 

rj+r'l+rl <4(r+-r_) 2 , (16) 

into the set of three already given in [17, sec. III]. The problem encountered, Eggeling indicated, stemmed from the 
fact that the third of the published constraints delimiting the domain of triseparability is of the third-degree. The 
added constraint (16) is necessary to ensure convexity. The plots in [17] other than Fig. 5 are themselves unaffected, 
he noted, since an appropriate numerical cutoff that took the "right" root of the third-degree polynomial was utilized 
in their generation. Fig. 1 is the same as Fig. 5 in [17], except that in addition to the central heart-shaped region 
of triseparable states displayed there, now there are also shown the three peripheral regions that the new fourth 
constraint (16) formally excludes. 



FIG. 1. Plot for the section r 




three peripheral regions that are excluded by the new fourth triseparability constraint (16). Without these three peripheral 
regions, the figure is identical to Fig. 5 in [9] 
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B. Probabilities for qubit case 



Based on some 654 million sampled points, we obtained in the three-qubit case, an estimate of .216769 that an EW- 
statc is biseparable (with respect to the partition 1|23) and .0630532 that it is triseparable. (As observed earlier, in this 
specific case, biseparability implies a positive partial transpose [17].) Table II is the qubit analogue of Table I. (Again, 
one billion random uniformly distributed points were examined in each randomly generated subsample to ascertain if 
they parameterized an EW-three-qubit density matrix. Approximately, 13.1 percent did in each subsample.) 



sample size 


Bures/SD prob. of biseparable state 


Bures/SD prob. of triseparable state 


130,866,541 


.216612 


.0631270 


130,860,088 


.216746 


.0630338 


130,871,093 


.216888 


.0631249 


130,852,647 


.216821 


.0630014 


130,867,972 


.216779 


.0629789 


654,326,341 


.216769 


.0630532 


std. dev. 


.0001026 


.000069 



TABLE II. Estimated Bures/SD probabilities, based on Monte Carlo simulations, that a three-qubit EW-state is: (1) 
biseparable with respect to the partition 1 1 23; and (2) triseparable. For each of the five subsamples, one billion points were 
initially generated. 
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So, the Bures/SD probabilities of biseparability and triseparability are considerably greater in the three-qubit case 
than in the threc-qutrit one. This is as would be expected on the basis of past analyses [1]. 



C. Supplementary Analyses 

Among the constraints that must be satisfied for an EW-state to be triseparable are [17, p. 6], 

0<r_<i; ^(1 - 2r_) < r + < 1 - 5r_. (17) 

If we integrate the normalized form of (15) over the full domain of parameters (1), except for imposing these constraints, 
we obtain an exact upper bound on the SD/Bures probability that a qutrit-qutrit-qutrit EW-state is triseparable of 

-70 - 325V2 esc" 1 V3 + 224^5 sin" 1 ,/f 

400^ " - 17766L < 18 > 

If we act similarly, but now set r_ = and use the normalized form of (14), we obtain an exact upper bound of 
27/64 = .421875 on the SD/Burcs probability that a qubit-qubit-qubit EW-state is triseparable. 

We can also get an apparently quite weak exact upper bound on the SD/Burcs probability that a three-qutrit 
EW-state is biseparable, by simply imposing the single constraint < r_ < 1/3 [17, p. 9]. This gives 

26V2 + 54csc-iV3 ^ ^ 

277T 

(while our Monte Carlo simulations [Table I] indicate that the actual Bures/SD probability is more on the order 
of .069). Though, somewhat to our disappointment, we have not been able to obtain exact Burcs probabilities 
of separability for the EW-states, T. Eggcling has indicated (App. II) that if one restricts one's consideration to 
the permutation- invariant EW states (which are commutative), then the required integrations can be performed 
analytically. 

As an additional investigative probe into the SD/Bures geometry of the EW-tripartite states, having U ® U ® U 
symmetry, we have estimated — first, for the three-qubit case — the ratio (« .34398) of the SD area of the boundary 
of the triseparable states to the SD area of the boundary of the biseparable states. To obtain this estimate, we 
extracted the 3x3 submatrix of the 5x5 metric tensor ((6)-(13)) corresponding to the variables n, r^ 1 and then 
took the square root (h) of its determinant (having, of course, set r_ =0, as is demanded in the d = 2 case). Then, 
we randomly generated values of ri,r2,r^ and computed for each of the various defining inequalities, in turn, that 
valuc(s) of r + which would saturate each. We then evaluated h if and only if that set of ri,r2,r^,r + also satisfied 
(without necessarily saturating) all the remaining required inequalities. The estimate of .34398 is based on some 
21,000,000 such points. (For the first 10,000,000 points, the estimate was .342709.) We have not explicitly ruled out 
the possibility that the SD areas themselves of the boundaries of the triseparable and biseparable states in question are 
unbounded in value, as preliminary analyses appear to indicate is the case for the area of the EW-states in general.) 

Similarly, for all the higher-dimensional thrcc-qudit cases (d > 2), where we now extract the 4x4 submatrix 
corresponding to r + , ri, r 2 , r 3 and solve for the values of r_ that saturate the various constraints, in turn, we obtained 
(on the basis of 18,000,000 points) an estimate of .0949602 for the ratio of the SD area of the boundary of triseparable 
states to the SD area of the boundary of biseparable states, and an estimate of .0000105263 for the ratio of the SD 
area of the boundary of triseparable states to the SD area of the boundary of states having positive partial transposes 
(with respect to the first tensor factor). (For the first 9,000,000 points, the estimates were, respectively, .0959525 and 
.000010494.) 



III. CONCLUDING REMARKS 



Let us indicate a number of papers of Batle et al [25-27]. In particular, these authors have investigated — following 
the review paper of Terhal [28] — the question of estimating the probabilities of (bi-) separability, using various 
separability criteria [27]. They have not employed the Bures measure in their various studies, as we have done here 
and elsewhere [1-4], but rather the (ZHSL) one utilized by Zyczkowski, Horodecki, Sanpera and Lewenstein in their 
pioneering paper [5]. However, we have observed in [29] that the ZHSL measure is not, in fact, proportional to 
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the volume element of any monotone metric, and, additionally, is based on an over-parametrization of the n 2 — 1- 
dimensional convex set ofnxn density matrices. So, one might argue, that the ZHSL measure serves best as a 
heuristic device, easily computable, but lacking a suitably rigorous theoretical rationale. 

Acin et al [18] have introduced a classification of mixed three-qubit states, in which they define the classes of 
separable, biseparable, Q and Greenberger-Horne-Zcilingcr states. Among other things, they conclude in contrast to 
the pure W-type states, that the mixed W class is not of measure zero since it is contained in a ball of finite radius (cf. 
[30,31]). Following the classification in [18], we have obtained in [4] — using quasi-Monte Carlo methods — estimates 
of lower bounds on the Burcs/SD probabilities of being a GHZ state and of being a W state. 

Clifton and Halverson [32] gave an elementary proof that the set of separable density matrices for any bipartite 
quantum system in which cither part is m/imfe-dimensional is trace-norm nowhere dense in the set of all density 
operators (cf. [33]). They assert that this result complements investigations [30] concluding that if both parts are 
finite-dimensional, then there is a (norm-dense) separable neighborhood of the maximally mixed state. However, in 
the case of 4 x 4 density matrices, if one adopts the maximal monotone metric, then it appears that the ratio of the 
volume of separable states to that of separable and nonseparable states is zero [4]. 
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IV. APPENDIX I OF T. EGGELING 



Computing the Bures distance of U <g) U ® [/-invariant states 



I start by recalling some known and useful facts about the nature of U <£> U <£> [/-invariant states. 

It is well known that the commutant of the set {U <g> U ® U} is spanned by the permutation operators of the 63 
More precisely every such symmetric state can be seen as a state (positive linear functional) on the finite dimcnsiona 
C*-algebra A consisting of the formal linear combinations of the group elements of G = 63 (the group algebr. 
C(G)). To associate to such a state u € A a density matrix p^ e B(C d ® C d ® C d ) we need a representatioi 
TTd'-Ai-^ B(C d <g> C d ® C d ). The most convenient is of course 



(1 



which extends easily to A. If we take for example d = 2 and 5 = (12) (3) we get the matrix representation 



V, 



(12)(3) 



ijk—1 



/10000000\ 
1 
1 
1 
1 
1 
1 
\ 1/ 

i.e. those doubly stochastic matrices acting on the kets of the canonical basis by swapping the tensor factors. A 
one may suggest this representation is not irreducible, i.e. it can be decomposed into irreducible subrepresentations 
Since A is spanned by 6 elements it is clear that the squares of the dimensions of the irreducible subrcprescntation 
(irrep) must sum up to 6. this is only possible as l 2 + l 2 + 2 2 or as l 2 + l 2 + l 2 + l 2 + l 2 + l 2 but since A is no 
abelian there must be at least one irrep of dimension greater than 1. This is exactly what lead us to the -R-operators 
In fact every element of A can be represented via two numbers corresponding to the two onedimcnsional component 
and one 2x2 matrix corresponding to the twodimcnsional irrep: 



p = r+ 



r 

n - 



-r 3 n - ir 2 
ir 2 r - r 3 



The full reduction of it is 



71 ~ n + ® n ~ ® ^0 



(2 



(3 



where 7r + is the trivial irrep assigning to each element always the number 1, 7r_ is the alternating irrep assigning t* 
each element in G its signum and n is the remaining twodimensional irrep. The ^s are the corresponding multiplicitic 
which are equal to the dimensions of the corresponding subspaces R + (C d (^)C d (^C d ) (Bose), i?_ (C d (^C d (^C d ) (Fermi 
and Ro(C d <g> C d £g> C d ) (para) divided by the dimension of the respective irrep (1 or 2). 

To compute the Bures distance we can now make use of the properties of such a representation. In fact a repre 
sentation is nothing but a *-homomorphism, i.e. a mapping respecting the *-operation (hermitian conjugation in ou 
case) and products: 



n(a*) = 7r(a)* n(ab) = Tr(a)n(b). 



For the Bures distance we get: 



d B (p,cr) 2 



2 - 2tr 



= 2-2tr[i?] + tr 



7> + 



= 2 - Iv^.frZsl- 2i>_*/rZs: 



for p = n d (r) and a = TV d (s), 



>rs\/r 



- 2tr [R]_ti 



\JVrsVr 



— 2vn tr \irn ( \l \frs\fr\ 1 



- 2tr [R] a tr 



>rs\/r 



(4 
(5 

rfi 



For the last summand we have to use the expression of Hiibner for 2x2 matrices: 



tr 



7TQ 



(r s + ns! + r 2 s 2 + r 3 s 3 ) + ^V r o - r i - r l ~ r l V s o - s i - s 2 - s 3- 



Taking r = 1 — r + — r_ into account we get for infinitesimal displacements: 

d B {p,p + dpf = 2- 2v + ^r + {r + +dr + ) - 2v_^r_{r_ + cZr_) 
1 



- 2f |- ((1 - r+ - r_)(l - r + - r_ - dr + - dr_) + n(ri + dri) + r 2 (r 2 + dr 2 ) + r 3 (r 3 + dr 3 )) + 
\-\J {I - r + - r_) 2 - r\ - rl - rl^f (1 - r + - r_ - dr + - dr_) 2 - {r x + dri) 2 - (r 2 + dr 2 ) 2 - (r 3 + dr 3 ) 2 



Differentiating ds{p,p + tdp) twice with respect to t at t = we get for the metric tensor: 

l(v_dr*_ v + dr\ ^ {dri + dr 2 + dr 2 + ^^^^T^f'^" 2 ) ' 



9ij(p)dr 7 dr 3 = 



4 \ r_ r + 



1 — r + — r_ 



For the volume clement we can thus derive 



^detg(p) = 



iy {v Q v + r- + v_(v Q r + - v + {\ - r + - r_))) 

2 ry* 2 /v>2 2 \ 



32(1 — r + — r_) y r + r_((l — r + — r_) 2 — r 2 — r 2 — r|) 



with 



v+ = 



d 3 + 3d + 2d 

~6 



d 3 - 3d 2 + 2d 
6 



VQ = 



d 3 -d 



Since for three Qubits (d=2) no completely antisymmetric states exist (the projector is trivial) we have to take v 
and r_ — for d=2. 



V. APPENDIX II OF T. EGGELING 



For commutative subsets of states (=abelian algebras) one can write down the Bures distance directly. Therefore, 
taking the subset of permutation invariant tripartite Werner states (which is commutative) we can compute the 
volume element and integrate over the triangles you see in the picture of my old PRA paper. This integration can be 
done analytically and leads to: 

p sep ,tr ip , P er m = -£:(-16 + 6^6 + 5 log 3(6 ~ ^ = 0.170502, (20) 

Pbise P ,tri P ,per m = ^(1 - 5^5 + 4V6 - 10 log | ± ) = 0.179607, (21) 

ppt=biseparable for permutation invariant tripartite states, so there is nothing to calculate. These numbers together 
with the Monte-Carlo simulations correspond to the intuition that triseparable and separable states are concentrated 
in the vicinity of the subset of permutation invariant states This is intuitive if one looks at the figures describing these 
sets of states. 



ACKNOWLEDGMENTS 



I would like to express appreciation to the Kavli Institute for Theoretical Physics for computational support in this 
research and to T. Eggcling for his sustained interest and involvement. 



[1] P. B. Slater, J. Phys. A 32, 5261 (1999). 

[2] P. B. Slater, Euro. Phys. J. B 17, 471 (2000). 

[3] P. B. Slater, Quant. Info. Proc. 1, 397 (2002). 

[4] P. B. Slater, quant-ph/0211150 (to appear in J. Phys. B). 

[5] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewenstein, 58, 883 (1998). 

[6] M. Hiibner, Phys. Lett. A 163, 239 (1992). 

[7] M. Hiibner, Phys. Lett. A 179, 226 (1993). 

[8] S. L. Braunstein and C. M. Caves, Phys. Rev. Lett. 72, 3439 (1994). 
[9] D. Petz and C. Sudar, J. Math. Phys, 37, 2662 (1996). 
[10] R. E. Kass, Statist. Sci. 4, 188 (1989). 

[11] B. R. Frieden, Physics from Fisher Information (Cambridge U. P., 1998). 
[12] H. P. Yuen and M. Lax, IEEE Trans. Infor. Th., 19, 740 (1973). 

[13] C. W. Helstrom, Quantum Detection and Estimation Theory (Academic, New York, 1976). 
[14] A. S. Holevo, Math. Methods Statist. 10, 6 (2001). 

[15] T. Tilma, M. S. Byrd and E. C. G. Sudarshan, J. Phys. A 35, 10445 (2002). 
[16] P. B. Slater, quant-ph/0308037. 

[17] T. Eggeling and R. Werner, Phys. Rev. A 63, 042111 (2001). 

[18] A. Aci'n, D. Bruss, M. Lewenstein and A. Sanpera, Phys. Rev. Lett. 87, 040401 (2001). 
[19] R. F. Werner, Phys. Rev. A 40, 4277 (1989). 
[20] J. Dittmann, J. Phys. A 32, 2663 (1999). 
[21] P. B. Slater, J. Geom. Phys. 39, 207 (2001). 
[22] R. Jozsa, J. Mod. Opt. 41, 2315 (1994). 

[23] J.-L. Chen, L. Fu, A. A. Ungar and X.-G. Zhao, Phys. Rev. A 65, 024303 (2002). 
[24] J. Dittmann, J. Geom. Phys. 31, 16 (1999). 

[25] J. Batle, A. R. Plastino, M. Casas and A. Plastino, J. Phys. A 35, 10311 (2002). 
[26] J. Batle, A. R. Plastino, M. Casas and A. Plastino, quant-ph/0306188. 
[27] J. Batle, A. R. Plastino, M. Casas and A. Plastino, quant-ph/0307167. 
[28] B. Terhal, Theor. Comp. Sci., 287, 313 (2002). 
[29] P. B. Slater, Lett. Math. Phys. 52, 343 (2000). 



10 



[30] S. L. Braunstein, C. M. Caves, R. Jozsa, N. Linden, S. Popescu and R. Schack, Phys. Rev. Lett. 83, 1054 (1999). 
[31] S. Szarek, quant-ph/0310061. 

[32] R. Clifton and H. Halverson, Phys. Rev. A 61, 012108 (2000). 
[33] P. B. Slater, J. Opt. B: Quantum Semiclass. Opt. 2, L19 (2000). 



11 



